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X-Ray binaries with neutron star accretors are traditionally divided in to two groups based 
on the masses of the donor stars. One is low-mass X-ray binaries (LMXBs) with donor stars 
less massive than ~ 1.5Mq, the other is high-mass X-ray binaries (HMXBs) with donor 
masses exceeding ^ 10. OM©. In LMXBs mass is exchanged through Roche-Lobe overflow 
(RLOF), while in HMXBs, the accretor is likely to be fed by the stellar wind-induced 
mass loss of the companion. Systems with donor masses between 1.5 and 10. DM© are 
called intermediate-mass X-ray binaries (IMXBs) . Few IMXBs have been discovered in the 
Galaxy. The reason is that, on one hand, mass transfer via RLOF is thought to be rapid 
and unstable due to the large mass ratio, leading to the formation of a common envelope; 
on the other hand, the donor stars at this mass range are unable to generate strong winds 



to power bright X-ray emission from the neutron star (jvan den Heuvel 



19751) . 



Recent investigations on IMXB evolution lead to important realization on the stability 
of super-Eddington mass transfer, and suggest that many, or perhaps most, of the current 
LMXBs descended from IMXBs. The studies on the evolution of the LMXB Cyg X-2 
Podsiadlowski fc RappaportjboOO : 



( King fc Ritter 



1998; 



Kolbetal 



I booo ) 



indicate that the 



mass of the donor star in this system must have been substantially larger (~ 3.5Mq) than 
its current value (~ 0.6Mq), implying that intermediate-mass systems can survive the 
high mass transfer phase by ejecting mo st of the transferred mass. The calculations by 



Tauris. van den Heuvel. fc Savoniid (j2000l ) have shown that the evolution of some IMXBs 
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may survive a spiral-in and experience a highly super-Eddington mass transfer phase on a 
(sub)thermal timescale if the convective envelope of the donor star is not too deep. These 
systems provide a new formation channel for binary millis econd pulsars with heav y CO 
white dwarfs and relatively short orbital periods (3—50 davs). lDavies fc HansenI (|l998r ) have 
independently suggested that IMXBs may be the progenitors of recycled pulsars in globular 
clusters. These works emphasize the necessity of accurately defining an evolutionary path of 
IMXBs, and motivate systematic analysis of binary systems undergoing thermal timescale 
mass transfer. 

Since the donor star in an IMXB is more massive than the accretor, the RL radius 
of the donor will shrink during the mass transfer. At the same time the donor star will 
either grow or shrink due to mass loss. The stabili ty of the mass transfer depends on the 
radiu s-mass exponents for the donor and its RL (jSoberman. Phinnev. fc van den Heuvel 



19971 ). ^2 = (aini?2/ainM2) and = (9 In i?L/a In M2)l, where 6 is the adiabatic or 



thermal response of the donor star to mass loss (M2 and R2 are the mass and radius of 
the donor star, Rl is its Roche- lobe radius, respectively). In general Rl decreases (^l > 0) 
when material is transfered from a relatively heavy donor to a light accretor, and vice 
versa. Donor stars with radiative envelopes will usually shrink (^2 > 0) in response to 
mass loss, while donor stars with a deep convective envelope expand rapidly (^2 < 0) in 
response to mass loss. The relative sizes of these parameters determine whether the mass 
transfer proceeds on either dynamical or thermal timescale. If ^2 > £,l the mass transfer is 
dynamically stable, occurring on either nuclear or thermal timescale. If ^2 < £,l the Roche 
lobe radius shrinks more rapidly than the adiabatic radius, and the mass transfer proceeds 



on dynamical timescale which leads to a common envelope and a spiral in phase. The final 
product cou ld be either a Thorn-Zytkow object or a short-period binary if the en velope is 
eje cted (e.g. 



Tauris. van den Heuvel. k, Savoniie 



Podsiadlowski et al 



200C; 



Podsiadlowski et al 



20021) . 



( 2002f ) made a survey of the X-ray binary sequences with the donor 



masses ranging from 0.6 to 7Mq. These authors found the actual mass transfer rates 
through RLOF sometimes deviate from the values given by the traditional formula for 
thermal mass transfeiQ, 

(iV/4 - M[) 



Mi 



th 



(1) 



where MJ and are the initial masses of the accretor and the donor respectively, and 
TKH is the Kelvin-Helmholtz time scale 



GMi 



(2) 



2R2L2 ' 

where G is the gravitational constant, and L2 the luminosity of the donor star. The work 
done by Langer et al. (2000) on the evolution of white dwarf binaries also indicates that 
Eq. (1) could overestimate the mean mass transfer rates by a factor of a few. 

However, due t o its simplicity, Eg . (1) has been widely u sed in population synthesis 



investigations (e.g. iHurlev et al 



2002 : 



Belczvnski et aLluOOT" ) for thermal timescale mass 



transfer in close binaries. The aim of this paper is to present a modified empirical formula to 
estimate the mean thermal timescale mass transfer rates onto neutron stars, by calculating 
the evolutions of IMXBs systematically. The results may be helpful to future investigations 
involving mass transfer processes in IMXBs. We describe the stellar evolution code and the 
binary model used in this study in §2. The calculated results and fitting formulae for the 
mass transfer rates are presented in §3. We conclude in §4. 



2. Binary calculations 



We have followed the evolution of the binary systems containing a neutron star and an 
intermed iate-mass second ary star using an updated version of the evolution code devel- 
oped by lEggletonl 



( 19711 



see also Pols et al. 1995). The opacit ies in the code are from 



Rogers fc IglesiasI (|l992l ) and from 



Alexander fc Ferguson 



(|l994l ) for temperatures below 
2^q3.8 assume a mixing length parameter of a = 2, and set the convective overshoot- 

ing parameter to be 0.2. The metallicity of the secondary is taken to be Z = 0.02 and 0.001, 
and the corresponding Helium abundance is 0.28 and 0.242, respectively. Each system is 
set to start with a neutron star of mass Mi — \AMq and a secondary of mass M2 from 1.6 
to 4.0Mq. Systems with donor mass higher than ^ 4.OM0 always experience unstable 



dyna mical mass transfer (jTauris. van den Heuvel. fc Savoniie 



2000; 



Podsiadlowski et al 



2002). The effe ctive RL radius of the secondary is calculated with Eggleton's equation 



(|Eggletonlll983l ). 



Rl 



0.49g' 



2/3 



0.6g2/3 + ln(l gi/3) ' 



(3) 



In 



Rapport. Di Stefano. fc SmithI ([199J) the thermal timescale mass transfer rate is taken to 



be half of that given by Eq. (1). 



where a is the orbital separation, and q — AI2/M1 is the mass ratio. We u se the following 
formula to calculate the mass transfer rate from the donor star via RLOF (|Eggletonlll97l[ ) 



-RMT{R2 - Rl?, 



(4) 



where RMT is a parameter adjusted automatically in the code, usually taken to be 500. 
The mass loss of the secondary via stellar w i nd is calculated according to the empirical 
formula given by iNieuwenhuiizen fc de Jageii (|1990[ ) , 



log(-M2w) = -14.02+ 1.241og(-^) + 0.16 log(—^; 

Lq IVIq 



0.811og(|^). 



(5) 



To follow the details of the mass transfer processes, we also include losses of orbital angular 
momentum due to mass loss, magnetic braking, and gravitational-wave radiation, although 
the last process is not important in this analysis. For magnetic braking, we use the stan- 
dard angular momentum prescription suggested bv lRapport. Verbunt. &: Josd (jl983r ). The 
Eddington luminosity of the neutron star is LEdd = 47rG'Mimp/tTT — 1-3 x 10^^(Mi/Mq) 
ergs~^, where nip and ctt are proton mass and the cross section of Thompson scattering, 
respectively. We limit the maximum accretion rate of the neutron star to the Eddington 
accretion rate MEdd — RiL-^dd/ GMi 1.5 x lO^^MQyr^^ for a 1.4Mq neutron star), 
i.e., Ml = — /M2, where 



/ 



(6) 



1 if -M2 < MEdd 

L -MEdd/M2 if -M2 > AfEdd ' 
We let the excess mass be lost from the system with the specific orbital angular momentum 
of th e neutron star. The orbital separation then ch anges according to the following equation 
(e.g. ISoberman. Phinnev. fc van den Heuvellll997l) 
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9(1-/), 



(7) 



2(1 + 9)^ 

where M — Mi + M2 is the total mass, J the orbital angular momentum, and Jmb the 
rate of orbital angular momentum loss by magnetic brakinig. 



3. Results 

We have calculated a large number evolutionary sequences for IMXBs with various initial 
donor mass and orbital period, so that mass transfer starts when the donor star is on 
early and late main sequence (cases al and a2), in the Hertzprung gap (cases bl, b2, b3) 
and on the giant branch (cases cl, c2 and c3), respectively. In Fig. 1 we show the initial 
distribution of the binaries in the AI2 vs. logo diagram. Triangles, squares, and diamonds 
in the figure correspond to the donor stars being main-sequence stars, subgiants, and red 
giants at the onset of mass transfer, respectively. We stop the calculation when either the 
mass transfer becomes dynamically unstable or dominated by the nuclear evolution of the 
donor. Obviously the occurrence of thermal timescale mass transfer depends on the initial 
mass ratio and the orbital period. When Z — 0.02 we find that steady case al to case 
b2 thermal timescale mass transfer is possible if the initial donor mass is between 1.6 and 
3.6M0, while systems either containing a more massive donor or in cases bS and c3 are 
subject to delayed unstable dynamical mass transfer. When Z = 0.001 the calculated results 



show that systems with donors between 1.8 and 3.6Mq experience stable thermal timescale 
mass transfer in cases al to c2, otherwise the mass transfer is dynamically unstable. 

Typical examples of the evolutionary sequences with Z = 0.02 and 0.001 are shown in 
Figs. 2 and 3 respectively. In Fig. 2 the system contains a neutron star and a companion 
with initial mass of S.OMq, which starts filling its RL roughly at the end of its central 
hydrogen burning. In Fig. 3 the donor has an initial mass of 3.6AfQ and starts filling its 
RL right after its helium ignition. In the figures we demonstrate the evolution of the mass 
transfer rate, the orbital period, the donor mass, and the neutron star mass with time. In 
Fig. 2 the mass transfer rate first rises rapidly to ~ 10~^'^ MQyr"^, then declines to a few 
10~* M0yr~^ after ~ 9 Myr, and stays around this value for ~ 1 Myr. During the former 
rapid mass transfer phase, the donor mass decreases from 3Mq to < IM©, but most of 
the mass is lost from the system, and efficient accretion by the neutron star occurs during 
the latter part of the mass transfer phase. The orbital period first decreases to around 1.2 
day, and then increases to ~ 30 day at the end of mass transfer. Mass transfer shown in 
Fig. 3 is more rapid due to the more massive and evolved donor star, lasting around 0.1 
Myr. About 2.6 Mq mass is transferred from the donor star during this phase, most of 
which is lost from the system, and the neutron star mass hardly changes. 

In our work the initiation (tj) and termination time (tf) of the thermal timescale mass 
transfer is assumed to be once the mass transfer rate exceeds and declines to the Eddington 
limit of the neutron star. The mean mass transfer rate Mmean is calculated from the fol- 
lowing equation, 

Mi, - Ml 

M^ean = (8) 
tf - ti 

where and Ml are the donor mass cAt = t\ and tf, respectively (stellar wind mass loss is 
negligible). In Tables 1 and 2 we list the calculated values of ti, tf, Mj, M|, Mmean and the 
maximum mass transfer rates M^ax. for evolutions with Z = 0.02 and 0.001, respectively. 
For comparison, we also list the expected values of Mth calculated with Eq. (1). Figures 4 
shows the calculated mean mass transfer rates as a function of (M2 — M{)/tkh- A linear 
fit can be obtained to be 

M_.0.28M^ (9) 
tkh 

for Z = 0.02, and 

(M| - Mi) 



Tkh 

for Z = 0.001. 



(10) 



4. Summary and discussion 

Our numerical calculations show that there are stable super-Eddington thermal timescale 
mass transfer processes in IMXB systems with donor mass between 1.6 and 3.6M0 in both 
Z = 0.02 and Z = 0.001 cases. We find that on average the traditional expression (Eq. [1]) 
overestimates the thermal timescale mass transfer ratcis by a factor of ^ 4. 

The results are obviously subject to various uncertainties in treating the mass trans- 
fer processes in binary evolution. One of the issues is the mass and angular momentum 



loss during mass transfer. We have used Eddinton-liniited accretion rate for neutron star 
accretion. Recent observations of quite a few binary millisecond radio pulsars constrain 



the pulsar masses to ~ 1.35 Mq (jBassa et al.l l20Q6l and references therein), suggesting 



that that almost all of the transferred mass may be lost rather accreted by the neutron 
star during the IMXB and LMXB phase. If the lost mass carries the specific orbital an- 
gular momentum of the neutron star, the orbital shrinking during thermal timescale mass 
transfer would be slower than we have calculated. This can be clearly seen from Eq. (7) 
by setting f — 0. The effect on mass transfer is most significant for those binaries with 
mass transfer rates being mildly supper-Eddington. For example, we find that the mean 
mass transfer rate is decreased by a factor ^ 5 if the donor mass is ~ 1.5 — 2.0 A/q and 
Z = 0.02. We note that in white dwarf binary evolution, it has also been realized that the 
strong mass loss from the accretor can stabilize the mass transfer even for a relatively high 



mass ratio, avoiding the formation of a commo n envelope (jHachisu. Kato. fc Nomoto 



Li fc van den Heuvel 



1997 



Langer et al 



circumbinary disk rather leaves the system (|Soberman. Phinnev. &: van den Heuvel 



1996 



20001). If, however, part of the lost mass forms a 



19971) 



the disk would extract orbital angular momentum from the binary by t idal torque, en 
hancing the mass transfer rates and leading to more rapid orbital shrink ij Spruit fc Taam 
20011) . 

It should be also noted that, during the evolution of an IMXB the strong X-ray ra- 
diation by the accretor and the accretion disk could illuminate the donor star and cause 



expansion of the donor and str ong stellar wind ijPodsiadlowski 



1991 



Hameury et al 



1993 



Phillips fc Podsiadlowskill2002r ). which would also lead to a higher mass transfer rate and 
shorter duration of mass transfer. As there is not a generally accepted theory on the irra- 
diation effect, we have not include it in our calculations. 
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Fig. 1. The initial distribution of donor mass and orbital separation at the onset of mass 
transfer. The triangles, squares and diamonds refer to donors on main sequence, subgiant 
branch and red giant branch, respectively. 
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Fig. 2. The evolution of an IMXB witli tlie donor of mass S.OMq and Z = 0.02. Tlie mass 
transfer starts riglit after the depletion of its core hydrogen. The solid and dotted lines show 
the evolution of the mass transfer rate and the donor mass in the upper panel, and the 
evolution of the orbital period and the neutron star mass in the lower panel, respectively. 
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Fig. 3. Similar as Figl. The evolution of an IMXB with the donor of mass 3.6Mq and 
Z = 0.001. The mass transfer starts right after its central helium ignition . The solid and 
dotted lines show the evolution of the mass transfer rate and the donor mass in the upper 
panel, and the evolution of the orbital period and the neutron star mass in the lower panel, 
respectively. 



10 




Fig. 4. The upper and lower panels show the fitting results for the mean thermal timescale 
mass transfer rates in Z = 0.02 and Z = 0.001 systems, respectively. The filled dots indicate 
the calculated results in this work. 
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Table 1. The calculated data of 24 evolutionary sequences with stable thermal timescale 
mass transfer {Z = 0.2). Mmean, -^max and Mth are the mean mass transfer rates, the max- 
imum mass transfer rates and the mass transfer rates calculated from Eq. (1), respectively. 



Ma 


Case 


ti 




tt 


Ml 


log Mmean 


log Mmax 


log Mth 


(Me) 




(10* yr) 


{Mq) 


(10* yr) 


{Mq) 


(Moyr-i) 


(Moyr-i) 


(Moyr-i) 


1.6 


al 


7.52 


1.58 


7.58 


1.51 


-7.94 


-7.90 


-7.72 


1.8 


al 


4.99 


1.80 


5.05 


1.63 


-7.61 


-7.55 


-7.27 


2.0 


al 


5.10 


2.00 


5.18 


1.61 


-7.33 


-7.14 


-6.87 


2.3 


al 


3.51 


2.30 


3.57 


1.52 


-6.91 


-6.54 


-6.52 


3.0 


al 


1.56 


3.00 


1.61 


1.26 


-6.45 


-5.70 


-5.99 


3.3 


al 


1.21 


3.30 


1.25 


1.13 


-6.23 


-5.48 


-5.81 


1.8 


a2 


11.0 


1.80 


11.1 


1.64 


-7.75 


-7.58 


-7.02 


2.0 


a2 


8.06 


2.00 


8.15 


1.60 


-7.35 


-6.97 


-6.70 


2.3 


a2 


5.31 


2.30 


5.39 


1.50 


-6.97 


-6.41 


-6.36 


2.6 


a2 


4.09 


2.61 


4.16 


1.36 


-6.73 


-5.99 


-6.04 


3.0 


a2 


2.68 


3.00 


2.73 


1.02 


-6.34 


-5.56 


-5.78 


3.3 


a2 


2.19 


3.30 


2.23 


0.95 


-6.17 


-5.35 


-5.56 


3.6 


a2 


1.73 


3.60 


1.76 


0.84 


-5.92 


-4.92 


-5.40 


1.6 


bl 


18.8 


1.57 


18.8 


1.52 


-7.99 


-7.96 


-7.31 


1.8 


bl 


13.2 


1.80 


13.2 


1.59 


-7.46 


-7.29 


-6.83 


2.0 


bl 


9.67 


2.00 


9.72 


1.52 


-7.02 


-6.68 


-6.44 


2.3 


bl 


6.46 


2.30 


6.51 


1.36 


-6.7 


-6.20 


-6.13 


2.6 


bl 


4.61 


2.60 


4.64 


1.03 


-6.34 


-5.83 


-5.84 


3.0 


bl 


3.14 


3.00 


3.16 


0.66 


-6.03 


-5.44 


-5.53 


3.3 


bl 


2.42 


3.30 


2.43 


0.65 


-5.75 


-5.22 


-5.39 


3.6 


bl 


1.92 


3.60 


1.92 


0.78 


-5.39 


-4.91 


-5.22 


3.0 


b2 


3.17 


3.00 


3.22 


0.40 


-6.30 


-5.44 


-5.46 


3.3 


b2 


2.47 


3.30 


2.49 


0.44 


-5.96 


-4.98 


-5.18 


3.6 


b2 


1.95 


3.60 


1.98 


0.46 


-6.00 


-4.96 


-5.10 
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Table 2. The calculated data of 34 evolutionary sequences with stable thermal timescale 
mass transfer {Z = 0.001). 



Ma 


Case 


ti 




ti 


Ml 


log Mmean 


log Mmax 


log Mth 


{Me) 




(10* yr) 




(10* yr) 


(Mo) 


(Moyr-^) 


(Moyr-i) 


(Moyr-i) 


1.8 


al 


3.29 


1.70 


3.47 


1.39 


-7.77 


-7.48 


-7.53 


1.8 


a2 


7.33 


1.80 


7.52 


1.21 


-7.50 


-6.23 


-7.05 


1.8 


b2 


9.01 


1.80 


9.13 


0.80 


-7.05 


-6.21 


-6.36 


2.0 


a2 


5.52 


2.00 


5.67 


1.17 


-7.27 


-6.40 


-6.89 


2.0 


bl 


6.40 


2.00 


6.55 


0.73 


-7.05 


-6.043 


-6.57 


2.0 


b2 


6.85 


1.87 


6.96 


0.69 


-6.99 


-4.14 


-6.08 


2.3 


cl 


6.04 


2.31 


6.05 


0.87 


-5.69 


-5.40 


-4.88 


2.3 


c2 


5.75 


2.31 


5.76 


0.96 


-5.82 


-5.64 


-5.03 


2.3 


c3 


6.10 


2.31 


6.11 


0.92 


-5.70 


-5.37 


-4.81 


2.6 


a2 


3.07 


2.62 


3.11 


0.84 


-6.36 


-5.50 


-6.13 


2.6 


bl 


3.35 


2.62 


3.38 


0.56 


-6.05 


-5.43 


-5.76 


2.6 


b2 


3.42 


2.62 


3.44 


0.47 


-6.01 


-5.31 


-5.58 


2.6 


cl 


4.44 


2.62 


4.442 


0.77 


-5.21 


-5.07 


-4.61 


2.6 


c2 


4.31 


2.62 


4.31 


0.73 


-5.20 


-4.91 


-4.76 


2.6 


c3 


4.43 


2.62 


4.43 


0.77 


-5.21 


-4.82 


-4.63 


2.8 


a2 


2.64 


2.79 


2.68 


0.81 


-6.26 


-5.39 


-5.83 


2.8 


bl 


2.93 


2.79 


2.99 


0.40 


-6.41 


-5.23 


-5.57 


2.8 


b2 


2.97 


2.79 


3.00 


0.42 


-6.07 


-5.10 


-5.28 


2.8 


cl 


3.80 


2.78 


3.80 


0.82 


-5.01 


-4.64 


-4.48 


2.8 


c2 


3.72 


2.79 


3.72 


0.79 


-5.18 


-4.71 


-4.61 


2.8 


c3 


3.76 


2.79 


3.76 


0.81 


-5.00 


-4.67 


-4.65 


3.0 


a2 


2.37 


3.03 


2.48 


0.40 


-6.63 


-5.13 


-5.59 


3.0 


bl 


2.40 


3.03 


2.47 


0.41 


-6.44 


-5.08 


-5.48 
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